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•'—J ■ Abstract 

^^' We calculate the open path phase in a two state model with a slowly (nearly 

adiabatically) varying time-periodic Hamiltonian and trace its continuous de- 
velopment during a period. We show that the topological (Berry) phase attains 
TT or 27r depending on whether there is or is not a degeneracy in the part of 
the parameter space enclosed by the trajectory. Oscillations are found in the 
*Cj . phase. As adiabaticity is approached, these become both more frequent and 

^S| ' less pronounced and the phase jump becomes increasingly more steep. Integral 

relations between the phase and the amplitude modulus (having the form of 
>0 ' Kramers-Kronig relations, but in the time domain) are used as an alternative 

^"l . way to calculate open path phases. These relations attest to the observable 

^— ^ ' nature of the open path phase. 

O 
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r-| ■ 1 Introduction 

In the last fifteen years much attention has been given to phases in wave functions 
and, in particular, to the topological (or Berry) phase, which is a signature of the 
trajectory of the system ([l]-[6]) and which is manifest in some interference and other 
5^ I experiments [7]. As noted in earlier works [1, 8, 9], the topological phase itvr that is 

picked up in a full revolution of the system is linked to the existence of a degeneracy 
of states (or crossing of potential energy surfaces) somewhere in the parameter space. 
This degeneracy need not be located in a region that is accessed to in the revolution; 
however, its removal even by a minute amount will cause the topological phase to be 
zero or an integral multiple of 27r. The physical model treated in this work confirms 
this effect; indeed the calculated topological phase shows a change from vr to 27r as 



the degeneracy disappears (cf. Figs 1 and 2). We tackle the problem by tracing a 
continuous variation in the non-cyclic, open path phase [4] (also named "connection" 
[2]), that is denoted in this work by 7(t) {t is time). 

To obtain an expression for 7(t) we study (for both the degenerate and non- 
degenerate alternatives) an explicitly solvable model. Both a detailed analysis and 
the figures exhibit, as a novelty, oscillations in 7(t). These become increasingly more 
frequent and of lesser magnitude as the adiabatic limit is approached. Furthermore 
it is observed that in the adiabatic limit of the degenerate case, the change in the 
open path phase is abrupt and results in a step function like behavior. 

In an alternative approach to the calculation of the open-path phase we develop 
reciprocal relations between phase and amplitude moduli of time dependent wave 
functions (Section 2). Versions of these relations in other contexts were given earlier 
([10]- [12]). The existence of these relations has the remarkable consequence that the 
associated open path phase, defined by them, is a "physical observable" (and inter 
alia gauge invariant) as a function of the path, a quality heretofore associated with 
the closed path (Berry) phase. 

2 Theory 

We start by invoking the Cauchy's integral formula which takes the form: 

where w{z) is analytic in the region surrounded by the anti-clockwise closed path. In 
what follows we choose the closed path to be the real axis t traversed in the reverse 
direction, of the infinite interval — oo < t < oo and an infinite semi circle in the lower 
half of the complex plane, as will be discussed later. We shall concentrate on the 
case that z is a real variable t and so equation (1) becomes: 

^ ^ vrz i_oo t' -t mJscC-t 

where P stands for the principal value of the integral, C = rexpiO, dC, = irexpiOdO 
and it is assumed that r ^ oo (the subscript SC in the second term stands for 
semi-circle). Next it is assumed that w{z) along the semi-circle is zero namely 

lim w{z) = 0, for 6' / 0, vr (3) 

Z — >OQ 

SO that equation (2) becomes: 

-Kl J-oo t' — t 

Assuming that the function w^z) is written as w^z) = w{t, y) = n(t, y) + iv{t, y) 
where z = t + iy, it can be shown by separating the real and the imaginary parts, 
that equation (4) yields the two equations: 

u{t) = --P r ^dt' and v{t) = -P r ^dt' (5) 

"^ J — OO T^ t TT J — oo ^ t 



These relations are of the Kramers-Kronig (KK) or dispersion equations type [13] 
and they will be applied in the time domain, u and v are Hilbert transforms [14]. 
Our aim is to employ equations (5) to form a relation between the phase factor in 
a wave function and its amplitude-modulus. If a wave- function amplitude ip{t) is 
written in the form: 

^(t)=f(t)exp(a(t)) (6) 

where T{t) and A(t) are real functions of a real variable t, the function w{z) which 
will be defined as: 

w{z)=ln{i){z)) = ln{t{z))+iX{z) (7) 

is assumed to fulfill the necessary conditions to employ the KK equations. This 
implies the following: (a) the function ^l^{z) is analytic and is free of zeroes in the 
lower complex half plane (however, ip{z) can have simple zeros on the real axis, as 
is made clear in several publications [10, 11, 14, 15]). (b) ip{z) becomes, along the 
corresponding infinite semi-circle, a constant (in fact this constant has to be equal 
to 1 but if the constant is 7^ 1 the analysis will be applied to ip{z) divided by this 
constant). Thus, identifying ln(r(t)) with u{t) and A(t) with v{t) we get from the 
second part in equation (5) the following expression: 

1 f^^ ln(f(t')) , , 

m = -P / ^r^dt' (8) 

Next assuming that r(f') is an even function, the equation for \{t) can be written 
as: 

and if T{t') is periodic then equation (9) can be further simplified to become: 



^<'>^?-rs-' 



\{t) = —P / dt'ln(rit')) y . (10) 

where T is the relevant period. 

3 The Model 

3.1 The Basic Equations 

In this work, the reciprocal relations in equation (5) are used in the form shown 
in equation (10) . A more general formulation of the reciprocal relations, includ- 
ing several applications, will be presented in a separate publication. Equation (10) 
is applied to two examples based on the Jahn- Teller model [16] which, following 
Longuet-Higgins [17], can be expressed in terms of an extended version of the Math- 
ieu equation, ([8, 18, 19]) namely: 

H = -\Eel^ - Gi{q, 4>) cos(2^) + ^2(9, 4>) M^O). (U) 



Here 9 is an angular (periodic) electronic coordinate, (j) is an angular nuclear periodic 
coordinate which is constrained by some external agent, as in ([6], [20]), to change 
linearly in time, namely cp = tot (thus, if T is the time-period, then to = i^), q is 
a radial coordinate, i?e« is a constant and Gi{q,(/)); i = 1,2 are two functions to be 
defined later. 

The Schrodinger equation is {h = 1): 

i_=«*. (12) 

and this will be solved approximately to the first order in ■^, for the case that the 
ground state is an electronic doublet. In a representation, adopted from [17], this 
doublet is described in terms of the electronic functions cos 9 and sin9 and therefore 
ip can be expressed as: ([6], [12]) 

^ = Xi{t)cos9 + X2{t)sm9 (13) 

In what follows equation (12) will be solved for the initial conditions: xi(t = 0) = 1 
and X2{t = 0) = 0. Replacing xi(0 ^^^ X2it) by ip+{t) and ip-{t) defined as: 

1 1 

V'±(i) = 2 exp(i-^e«t)(xi T iX2) (14) 

we get the corresponding equations for ip+{t) and %lj^{t): 

itpj^ = — GV'- and iijj- = — G*'(/'+ (15) 

where G is defined as: G = Gi + iG2, and the dot represents the time derivative. 

Next we eliminate ip- from equation (15) to obtain a single, second order equa- 
tion for ip^: 

V;+-ln(G)V'+ + ^|G'|V+ = 0. (16) 

Writing G = |G|exp(i<^) we shall be interested in cases where |G| is constant, so 
that only <I> is time-dependent. Thus equation (16) becomes: 

V^+-i6^+ + ^|G|V = 0. (17) 

Once equation (17) is solved we can obtain xi(^)) the eigen-function for the initially 
populated state. Usually, this is a fast oscillating function of t where the oscillations 
are caused by the "dynamical phase" ^IGjt. This oscillatory component is eliminated 
upon multiplying Xi(*) by ex-p{—^i\G\t). In what follows we consider the smoother 
function r]{t) defined as: 

7?(t)=Xi(t)exp(-ii|G'|t) (18) 

Our aim is the study of the time dependence of the phase 7(f) defined through the 
expression: 

r?(t) = p(f)exp(i7(t)) (19) 



with p{t) and 7(t) real. Once rj{t) is derived there are several ways to extract 7(t) 
we shall use the following two: (a) the first is the following analytical representation 
of the open path phase given by Pati [4] : 

7(t) = 9(ln(r/(t))) (20) 

where 9 stands for the imaginary part of the expression in the parentheses. Equation 
(20) is used for analytical purposes, as presented below. Special emphasis will be 
put on 7(t) at t = T where T is the period of the external field. The case of an 
arbitrary T will be discussed only briefly and we will be mainly interested in the 
adiabatic case where T is large, namely T S> \G\'^ for which ^(t = T) becomes the 
topological (Berry) phase (3. In what follows we distinguish between two cases: (a) 
The degenerate case for which the functions in equation (11) are given by: 

Gi((?,0) = Gqcos{(t>), G2{qA) = Gqsm{(t>). (21) 

Where G is constant. We term it the degenerate case because the two lowest eigen- 
values of equation (11) become equal in the (g, (f)) plane at (7 = 0. It is also noticed 
that: G = Gq. (b) The non-degenerate case. This is characterized by the condition 
that Gi = and G2 = cannot be simultaneously satisfied for real q and </>.It is 
not trivial to achieve this by a simple change of the expressions in equation (21) , 
since, e.g., adding a constant will only displace the real root, as has been previously 
discussed [19]. However, non-degeneracy can be attained upon replacing G2 by a 
quadratic polynomial in qsin^cp), such that the polynomial has no real roots. A term 
of this form is physically realizable in a low-symmetric molecular environment. (In a 
realistic case of non crossing potential energy surfaces for NaFH, the expression con- 
structed for G2 is very complicated [21]). Unfortunately, the equation (16) cannot 
be solved analytically for a general polynomial G2. Below (in section 4) we present 
an approximate solution for a case that a degeneracy is encountered neither at q = 
nor at any other real g-value. This is achieved by the choice: 



Giiq,^) = Gqcosi^), and G2(g, 0) = V(Gg)2 sin^ (,/.)+ ^2 (22) 

The quantity fi is related to the separation between the two potential energy surfaces. 
If we expand the square root for small q, we indeed recover a quadratic polynomial 



approximation, that has no real roots. It is noticed that now |G| = yGf + G2 = 

\J{GqY + M^- In what follows it is assumed for simplicity that the particle trajectory 
is on the circle 5=1. 

3.2 The Degenerate Case 

We start by considering the degenerate case and therefore in equation (16) $ = (/> = 
u)t and |G| = G as already mentioned, (see equation (21) ). As a result, equation 
(17) becomes: 

^+ - iuji:+ + ^GV+ = 0. (23) 



The solution of this equation (as well as that of a similar equation for tp^{t)) can be 
written in terms of trigonometric functions. Returning to the original %- functions 
we get for Xi{t) the following explicit expression: 

Xi = cos{kt) cos{— cot) -\ sm(kt) sm{—ujt) 

Li ZiK Zi 

G 1 

+ i— sin(/ct)cos(-a;t) (24) 



where /c, defined as: 



1 



fc = -v/G2+w2 (25) 

forms, together with w, two characteristic periodicities of the system. 

In Fig. 1 are shown several 7(t)-functions as calculated for three different values 
of G and T{= -^). It is noticed that as T increases, namely as the adiabatic limit is 
approached, ^{t) tends to a step function and (5- the Berry phase- reaches the value 
of vr. This behavior is also derived analytically as follows. 

Considering the case that T ^ oo (or uj ^ 0), one can show employing equations 
(20) and (24) for the adiabatic case, that ^(t) takes the form (discarding second 
order terms in w): 

hm (7(t)) = 9{ln[cos(i^t) + 0(a;)]} (26) 

I — ^oo Z 

Having this expression it is recognized that since cos(ia;t) > for t < -^ and 
cos{^ujt) < for t > 1^ it follows that -/{t) ~ for < t < |^ and j{t) ~ vr for 
^ < t < T. This also implies that the topological (Berry) phase /? ~ vr. From Fig. 
1 it is noticed that, when the adiabatic limit is approached (namely, T -^ oo), 7(t) 
becomes a step function. The step takes place at t ~ -g. It is therefore of interest 
to study the behavior of j{t) in the vicinity of t = 2"- Thus expanding j{t) around 
this value and keeping only first order terms in (t — -g) yield: 

T T I 

7(t Rs — ) = 9{ln[ t+ -smkt exp(jfct)]} (27) 

It is noticed that around t = n" ^^^ phase factor 7(t) oscillates (Rabi oscillations) 
and its periodicity is (-^). These oscillations become more frequent the larger is the 
value of the product GT(^> 1). 

In order to obtain the phase using equation (10) we have to construct from 
Xi(i), which when analytically continued to the complex plane becomes Xi{^)^ ^ 
new function that fulfills the requirements imposed on ip{z). The complex function 
Xi{z) is obtained by replacing in equation (24) , the variable t by z defined as: 

z = rexp(i(9); where < 6* < 27r and r > 0. (28) 

The first requirement imposed on 1^(2) is that it does not have zeros in the lower 
half plane. The newly formed function Xi(-2) has, in general, zeros in the lower half 
plane. But we have been able to show generally that near the adiabatic limit there are 
no zeros for the ground state in the lower half plane. [Moreover, a detailed numerical 
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Figure 1: The time-dependent phase factor 7 as a function of time calculated for 
the degenerate case. The calculations were done for different values of the external 
field frequency u}{= ^), and the coupling intensity G. In all three cases the Rabi 
oscillations are of a time period (-^) where k is defined as A; = 2 VC^ + w^. In each 
sub-figure two curves are shown; one, drawn as a full line, is the curve calculated 
employing equation (24) and the other, drawn as a dashed line, is the approximate 
curve calculated employing the Kramers-Kronig relation shown in equations (8-10) 
and (32). (a) T = 10^, G = 0.01. (b) T = 10^, G = 0.02. (c) T = 10^ G = 0.01. 



study showed that when the ratio of inverse periods (— ) = integer, the zeros (of the 
ground state) are located in the upper half plane (including the real axis). For the 
near adiabatic situation where (— ) is large, the requirement — = (a large) integer can 
(on physical grounds) differ only insignificantly from neighboring values of — that are 
non integral. We thus have two independent reasons for the assertion regarding the 
location of zeros in the near adiabatic case. This is also confirmed by our graphical 
results in Figures 1 and 2, which clearly show the increasing validity of the integral 
relations, as the adiabatic limit is approached, upon going from (a) to (c), and this 
even though the ratio (— ) is not chosen to be an integer.] The second requirement 
imposed on ip^z) is that it becomes equal to 1 along the infinite semi-circle on the 
lower half of the complex plane. From equation (24) it is readily seen that for 
T ^ oo the function Xi(-2) in the adiabatic limit becomes (for 9 > tt): 

lim XI (z) = I exp(i(fc + loj)Te'^) (29) 

T— >oo 2 2 

or 

lim xi(-2) = — exp(— (/c H — uj)t sin 9) exp{i(k -\ — uj)tcos9) (30) 

r->oo 2 2 2 

Therefore multiplying Xi(-2) by 2exp(— i(/c+ ^u;)re*^) yields the function ip{z) which 
becomes equal to 1 along the infinite-semi circle. Thus the function to be employed 
in equation (10) is ■ip{t) defined as: 

i^it)=2xiit)expi-i{k+^u;)t) (31) 

Combining Eqs. equation (8) , equation (10) , equation (18) , equation (19) and 
equation (31) we obtain the final expression for the phase 7(t) up to a linear function 
of time that follows from the KK equations: 

j{t) = —P / dt'ln(r{t')) y = (32) 

where T{t) is is absolute value of Xi{t)- It is important to emphasize that T is 
not necessarily equal to T (in our particular case T is equal to 2T). In Fig. 1 is 
presented 7(t) also as calculated from equation (32) . The results along the interval 
< i < (-2) were taken as they are but those along the interval -g- < i < 7" were 
found to be 27r below the values obtained by the direct method. We added to each 
of the calculated values the physically unimportant magnitude 27r. The comparison 
between the curves due to the two different calculations reveals a reasonable fit 
which improves when either T or G become large enough, namely upon approaching 
the adiabatic limit. Even the (Rabi) oscillations at the near adiabatic limit are well 
reproduced by the present theory. Moreover the theory yields the correct geometrical 
phase. It is also important to mention that when we are far from the adiabatic limit 
the fit is less satisfactory. However, we also found that for the choices of k which make 
the function Xi(^) periodic, namely when (— ) = integer, the agreement resurfaces 
[12]. 



4 The Non-Degenerate Case 

This arises when ^ 7^ (see equation (22) ). As a result we obtain for jGI, <& and $ 
the following expressions: 



sin Lut 



\G\ = JG^ + fi"^; ^ = arccos{p cos ujt); ^ = puj—^=^==== (33) 

^ V 1 — p^ cos^ ujt 

where p is defined as p = , In what follows we consider only the case when 

H is small enough so that |G| and <I> are, as before, equal to G and cot, respectively, 
but <i> will be written as: 6 = cj i^'."^*! . Thus equation (17) becomes: 

V^+tW+ + ^|G'|V+ = 0. (34) 

where the minus sign is for the < t < ( 2") - the first half period and the plus sign 
for -2 ^ t < T - the second half. For the first half we have the same equation as 
before and therefore also the same solution (see equation (24) ). As for the second 
half period we obtain a somewhat more complicated expression for the solution due 
to the matching of the two solutions at t = -E. Thus: 

Xi{t) = e"*''{cos(A;t) cos{-Lot) - -^ sm{k{T - t)) sin(-wt) 

G 1 

+ i—siTi.{kt)cos{-ujt)} 

In order to obtain the phase factor for the adiabatic case, equation (20) is applied 
as before, where Xi(*) is given by equation (35) . We employed equation (35) to 
calculate 7(t) based on the KK dispersion relations. 

In Fig. 2 are presented the results due to the two types of calculations as obtained 
for three sets of values of the parameters G and T. It is noticed that when either 
T or G become large enough (namely, approaching the adiabatic limit), as in the 
previous case, a reasonably good fit is obtained between the results due to the direct 
calculations and the ones based on the KK relations (equations (10) and (35)). 
Moreover, in this case, too, this new formalism yields the correct geometrical phase. 

The same analytic treatment can be done for the non-degenerate two-state model. 
Considering again the case that T ^ 00 (or a; ^ 0), but for equation (35) , we obtain 
that 7(t) takes the form: 

hm (7(t)) = 9{ln[cos(iu;t) + 0{u)]} + e{t - ^)7r (36) 

where Q{x) is the Heavyside function defined as being equal to zero for x < and 
equal to 1 for a; > 0. It is noticed that the sign of the expression in the square 
brackets is positive for < t < ^ which means that the phase factor is altogether 
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Figure 2: The time dependent phase factor 7 as a function of time calculated for 
the non-degenerate case. Details as in Fig.l. The approximate curve was calculated 
employing equations (8-10) and (32) for < t < (-g) and equations (8-10) and (35) 
for I < t < T. 
G = 0.01. 



(a) T = 10^, G = 0.01. (b) T = 10^ G = 0.02. (c) T = 10^, 
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zero (because also 0(x) = 0) but the sign is positive for ^ < t < T and therefore 
altogether 7(i) = 27r and this leads to a topological (Berry) angle P = 2tt. This 
result is expected because the Berry phase has to be zero (or 27r) in the case of no 
degeneracy. 

5 Conclusions 

The expression of the degeneracy and near-degeneracy dichotomy in the topological 
phase is the main subject of this paper. The respective values of vr and 2it after 
one revolution (seen in Figures 1 and 2, respectively) obtained in a two-stage model 
confirm the expectations. However, on the way to this result we earned some new 
results and insights. Oscillations near the half period stage were found (equation 
(27) ) and explained. We also studied the tendency of this and of other features 
in the "connection" (namely the non-cyclic phase) with the approach to adiabatic 
(slow) behavior. 

An attempt has been made in this article to establish a link between the time 
dependent phase (and its particular value, the topological phase, after a full rev- 
olution) with the corresponding amplitude modulus. To establish this relation we 
considered two alternative two-state models, exposed to an external field, under 
adiabatic and quasi-adiabatic conditions ([3], [19]). The two types of models are 
physically different: (a) one model contains an (ordinary Jahn- Teller type) degen- 
eracy at a point in configuration space; (b) the second is characterized by a nearly 
(in fact, non-) degenerate situation (of the pseudo-Jahn- Teller type([16])) where the 
two eigenvalues approach each other at some point in configuration space but do 
not touch. In Figs. 1 and 2 are presented time dependent phases and the (Berry) 
topological phases for these two models calculated in two different ways: once di- 
rectly by employing equation (20) and once by using the KK relations which led to 
equation (10) . Essentially these findings suggest that one may be able to obtain 
the time dependence of the phase from a series of time-dependent measurements of 
relative populations of a given state. We end by offering the following interpretation 
for our findings: The phase on the left-hand side of equations (8-10) is not a "phys- 
ical observable" in the conventional sense since no hermitian operator is associated 
with it [10]. Yet, phases have been observed in interference and other experiments 
[7]. In the present formulation, equations . (8-10) associate the observable phase of 
the wave function with T(t) (the observable probability amplitude) through integral 
expressions, in a similar way to that done in Ref. [11] for radiation fields. 
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